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Abstract. Some of recent important developments are overviewed, several 
longstanding open problems are discussed, and a perspective is presented for 
the mathematical theory of multidimensional conservation laws. Some basic 
features and phenomena of multidimensional hyperbolic conservation laws are 
revealed, and some samples of multidimensional systems/models and related 
important problems are presented and analyzed with emphasis on the proto- 
types that have been solved or may be expected to be solved rigorously at 
least for some cases. In particular, multidimensional steady supersonic prob- 
lems and transonic problems, shock reflection-diffraction problems, and related 
effective nonlinear approaches are analyzed. A theory of divergence-measure 
vector fields and related analytical frameworks for the analysis of entropy so- 
lutions arc discussed. 
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1. Introduction 

We overview some of recent important developments, discuss several longstand- 
ing open problems, and present a perspective for the mathematical theory of mul- 
tidimensional (M-D, for short) conservation laws. 

Hyperbolic Conservation Laws, quasilinear hyperbolic systems in divergence form, 
are one of the most important classes of nonlinear partial differential equations 
(PDEs for short). They typically take the form: 

d t u + V x • f (u) = 0, uel m , (l.i) 

for (i,x) e R+ +1 := R+ x R d := [0,oo) x R d , where V x = (d Xl , . . . , d Xd ), and 
f = (fi, . . . ,f d ) : K m -> (R m ) d is a nonlinear mapping with i t : R m -> R m for 
i = 1, . . . , d. Another prototypical form is 

d t A {u t , V x u) + V x • A(u*, V x u) = 0, (1.2) 

where A {p ,p) : RxR d -> M and A(p ,p) : MxM d -> R d with p = ( Pl ,...,p d ) e R d 
are C 1 nonlinear mappings. 

The hyperbolicity for (1.1) requires that, for any n £ s d -\ 

(V w f (w) • n) mxm have m real eigenvalues Aj(w; n), 1 < i < m. (1.3) 

We say that system (1.1) is hyperbolic in a state domain V if condition (1.3) holds 
for any wdP and n £ iS d_1 . 

The hyperbolicity for (1.2) requires that, rewriting (1.2) into second-order non- 
divergence form such that the terms u XiXj have corresponding coefficients a^, i, j = 
0, 1, d, all the eigenvalues of the matrix (aij)dxd are real, and one eigenvalue has 
a different sign from the other d eigenvalues. 

An archetype of nonlinear hyperbolic systems of conservation laws is the Euler 
equations for compressible fluids in R d , which are a system of d + 2 conservation 
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laws: 

dtp + V x • m = 0, 

d t m + V x -(^)+V x p = 0, (L4) 
dt(pE) + V x -(in(E + p/p))=0. 
System (1.4) is closed by the constitutive relations: 

1 Iml 2 

p = p(p,e), E=- l —L + e. (1.5) 

2 p A 

In (1.4)-(1.5), r = 1/p is the deformation gradient (the specific volume for fluids, 
the strain for solids), m = (mi, . . . , md) T is the fluid momentum vector with v — 
m/p the fluid velocity, p is the scalar pressure, and E is the total energy with e 
the internal energy, which is a given function of (r,p) or (p,p) defined through 
thermodynamical relations. The notation a ® b denotes the tensor product of the 
vectors a and b. The other two thermodynamic variables are the temperature 9 and 
the entropy S. If (p, S) arc chosen as independent variables, then the constitutive 
relations can be written as (e,p,9) — (e(p, S),p(p, S),9(p, S)) governed by the 
Gibbs relation: 

9dS = de + pdr = de -dp. (1.6) 

For a poly tropic gas, 

7-1 -£ 

wherc R > may be taken as the universal gas constant divided by the effective 
molecular weight of the particular gas, c v > is the specific heat at constant 
volume, 7 > 1 is the adiabatic exponent, and k > can be any positive constant 
by scaling. The sonic speed is c := ^p p (p, S). For polytropic gas, c = \fypjp- 

As indicated in §2.4 below, no matter how smooth the initial function is for the 
Cauchy problem, the solution of (1.4) generically develops singularities in a finite 
time. Then system (1.4) is complemented by the Clausius-Duhem inequality: 

dt(pS) + V x • (m5) > (1.8) 

in the sense of distributions in order to single out physical discontinuous solutions, 
so-called entropy solutions. 

When a flow is isentropic (i.e., the entropy S is a uniform constant So in the 
flow), the Euler equations take the following simpler form: 

( d t p + V x • m = 0, 

(1.9) 



p = Rp9 = np^e^, e = c v 9 = — — p 7 " 1 ^, 7 = 1 + — , (1.7) 

7-1 C„ 



9 t m + V x - (=fsi) +V xP = 0, 



where the pressure is regarded as a function of density, p = p(p, So), with constant 
So- For a polytropic gas, 

p{p)=K P \ 7 >1, (1.10) 
where n > can be any positive constant under scaling. When the entropy is 
initially a uniform constant, a smooth solution of (1.4) satisfies the equations in 
(1.9). Furthermore, it should be observed that the solutions of system (1.9) are 
also close to the solutions of system (1.4) even after shocks form, since the entropy 
increases across a shock to third-order in wave strength for solutions of (1.4), while 
in (1.9) the entropy S is constant. Moreover, system (1.9) is an excellent model for 
the isothermal fluid flow with 7 = 1 and for the shallow water flow with 7 = 2. 
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An important example for (1.2) is the potential flow equation: 

dtp^t, V x $) + V x • (p(* t) V X $)V X $) = 0, (1.11) 

where p($ t , V x $) = /i -1 ^ - (9 t $ + ±|V X <I>| 2 )), ft'(p) = p'(p)/p, for the pressure 
function p(p), which typically takes form (1.10) with 7 > 1. 

The importance of the potential flow equation (1.11) in the time-dependent Eulcr 
flows (1.4) with weak discontinuities was observed by Hadamard [123]. Also see Bers 
[16], Cole-Cook [81], Courant-Friedrichs [84], Majda- Thomas [159], and Morwawetz 
[165]. 

2. Basic Features and Phenomena 

We first reveal some basic features and phenomena of M-D hyperbolic conserva- 
tion laws with state domain T>, i.e., ueD. 

2.1. Convex Entropy and Symmetrization. A function i] : T> — > M is called 
an entropy of system (1.1) if there exists a vector function q : T> — > M d ,q = 
(qi,...,q d ), satisfying 

Vq i (u)=V/ l (u)Vf i (u), i = l,...,d. (2.1) 

Then q is called the corresponding entropy flux, and (77, q) is simply called an 
entropy pair. An entropy 77(11) is called a convex entropy in V if V 2 77(u) > for 
any u £ T> and a strictly convex entropy in V if V 2 r;(u) > cqI with a constant 
Co > uniform for u € X>i for any 2?i c 2?i <s 2?, where 7 is the m x m identity 
matrix. Then the correspondence of the Clausius-Duhem inequality (1.8) in the 
context of hyperbolic conservation laws is the Lax entropy inequality: 

d tV (u) + V x • q(u) < (2.2) 

in the sense of distributions for any C 2 convex entropy pair (77, q). 

The following important observation is due to Friedrich-Lax [109], Godunov 
[117], and Boillat [21] (also see Ruggeri-Strumia [179]). 

Theorem 2.1. A system in (1.1) endowed with a strictly convex entropy inT> must 
be symmetrizable and hence hyperbolic in T>. 

This theorem is particularly useful for determining whether a large physical 
system is symmetrizable and hence hyperbolic, since most of physical systems from 
continuum physics are endowed with a strictly convex entropy. In particular, for 
system (1.4), 

(77*, q*) = (-pS,-mS) (2.3) 

is a strictly convex entropy pair when p > and p > 0; while, for system (1.9), the 
mechanical energy and energy flux pair: 

, ,1 Iml 2 N 1 Iml 2 , , pip),, 

^*,q* = p (^—r + e > m b4" + e W + — 2.4 
2 p z 2 p z p 

is a strictly convex entropy pair when p > for polytropic gases. For a hyperbolic 
system of conservation laws without a strictly convex entropy, it is possible to 
enlarge the system so that the enlarged system is endowed with a globally defined, 
strictly convex entropy. See Brenier [24], Dafermos [91], Demoulini-Stuart-Tzavaras 
[94], Qin [176], and Serre [185]. 
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There are several direct, important applications of Theorem 2.1 based on the 
symmetry property of system (1.1) endowed with a strictly convex entropy. We list 
three of them below. 

Local Existence of Classical Solutions. Consider the Cauchy problem for a 
general hyperbolic system (1.1) with a strictly convex entropy: 



Theorem 2.2. Assume that u : R d -> V is in H s n L°° with s > d/2 + 1. 
Then, for the Cauchy problem (1.1) and (2.5), there exists a finite time T = 
T(||uo|| s , ||uo||l°°) € (0, oo) such that there exists a unique bounded classical so- 
lution u e C a ([0,T] x R d ) with u(i,x) e 2? /or (i,x) e [0,T] x R d and u <= 
CdCTjj^nC^ICT];^- 1 ). 

Kato [131] first formulated and applied a basic idea in the semigroup theory to 
yield the local existence of smooth solutions to (1.1). Majda in [157] provided a 
proof which relies solely on the elementary linear existence theory for symmetric 
hyperbolic systems with smooth coefficients via a classical iteration scheme (cf. [85]) 
by using the symmetry of system (1.1). Moreover, a sharp continuation principle 
was also provided: For u € H s with s > d/2 + 1, the interval [0,T) with T < oo 
is the maximal interval of the classical H s existence for (1.1) if and only if cither 
u(t, x) escapes every compact subset K <e V as t — > T, or ||(u t , Du)(t, OIU 00 — > oo 
as t — > T. The first catastrophe is associated with a blow-up phenomenon such as 
focusing and concentration, and the second in this principle is associated with the 
formation of shocks and vorticity waves, among others, in the smooth solutions. 
Also see Makino-Ukai-Kawashima [160] and Chemin [37] for the local existence of 
classical solutions of the Cauchy problem for the M-D Eulcr equations. 
2.1.2. Stability of Lipschitz Solutions, Rarefaction Waves, and Vacuum 
States in the Class of Entropy Solutions in L°° . Assume that system (1.1) 
is endowed with a strictly convex entropy on compact subsets of V. 

Theorem 2.3 (Dafermos [90, 91]). Suppose that w is a Lipschitz solution of (1.1) 
on [0,T), taking values in a convex compact subset K ofV, with initial data w . 
Let u be any entropy solution of (1.1) on [0,T), taking values in K, with initial 
data u . Then 



holds for any R > and t S [0,T) ; with L > depending solely on K and the 
Lipschitz constant o/w. 

The results can be extended to M-D hyperbolic systems of conservation laws with 
partially convex entropies and involutions; see Dafermos [91] (also see [22]). Some 
further ideas have been developed [42] to show the stability of planar rarefaction 
waves and vacuum states in the class of entropy solutions in L°° for the M-D Euler 
equations by using the theory of divergence- measure fields (see §7). 

Theorem 2.4 (Chen-Chen [42]). Let n e S^ 1 . Let R(i,x) = (/5,m)( x fi) be a 
planar solution, consisting of planar rarefaction waves and possible vacuum states, 
of the Riemann problem: 

R|t=o = (p±,m±), ±x-n>0, 



u| t=0 = u . 



(2.5) 
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with constant states (p±,m±). Suppose that u(t,x) = (p, m)(i,x) is an entropy 
solution in L°° of (1.9) i/iat moj/ contain vacuum. Then, for any R > and 
i € [0,oo), 

/ a(u,R)(i,x)dx < / a(u,R)(0,x)dx, 

J|x|<K J\x\<R+Lt 

where L > depends solely on the bounds of the solutions u and R, and a(u, R) = 
(u — R) T (/ Q L V 2 77*(R + t(u — R))dr)(u — R) with the mechanical energy n*(u) in 
(2.4). ' 

A similar theorem to Theorem 2.4 was also established for the adiabatic Euler 
equations (1.4) with appropriate chosen entropy function in Chen-Chen [42]; also 
cf. Chen [41] and Chcn-Frid-Li [55]. Also see Perthamc [174] for the time-decay 
of the internal energy and the density for entropy solutions to (1.4) for polytropic 
gases (also cf. [37]). 

2.1.3. Local Existence of Shock-Front Solutions. Shock-front solutions, the 
simplest type of discontinuous solutions, are the most important discontinuous non- 
linear progressing wave solutions of conservation laws (1.1). Shock- front solutions 
are discontinuous piecewise smooth entropy solutions with the following structure: 

(i) There exist a C 2 space-time hypersurface S(t) defined in (i, x) for < t < 
T with space-time normal (n t ,n x ) = (n t ,ni, . . . , n<j) and two C 1 vector- 
valued functions: u ± (t, x), defined on respective domains V ± on either side 
of the hypersurface S(t), and satisfying 

9tU ± + V-f(u ± ) = inV ± ; (2.6) 

(ii) The jump across S(t) satisfies the Rankine-Hugoniot condition: 

(n t (u+ - u") + n x • (f (u+) - f (u-))) \ s = 0. (2.7) 

Since (1.1) is nonlinear, the surface S is not known in advance and must be deter- 
mined as a part of the solution of the problem; thus the equations in (2.6)-(2.7) 
describe a M-D free-boundary value problem for (1.1). 

The initial functions yielding shock-front solutions are defined as follows. Let 
So be a smooth hypersurface parameterized by a, and let n(a) = (ni, . . . , n^)(a) 
be a unit normal to <So- Define the piecewise smooth initial functions u (x) for 
respective domains on either side of the hypersurface So- It is assumed that the 
initial jump satisfies the Rankine-Hugoniot condition, i.e., there is a smooth scalar 
function a (a) so that 

-a(a)(u+(a) - Uq (a)) + n(a) • (f(u+(a)) - f(u^(a))) = (2.8) 

and that a(a) does not define a characteristic direction, i.e., 

a(a) ^ Aj(u^), ae~S~o~, 1 < i < m, (2.9) 

where A,, i = 1, . . . , m, are the eigenvalues of (1.1). It is natural to require that 
S(0) = So. 

Consider the 3-D full Euler equations in (1.4) for u = (p, m, pE), away from 
vacuum, with piecewise smooth initial data: 

u| t=0 = u±(x), xeZtf. (2.10) 
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Theorem 2.5 (Majda [156]). Assume that So is a smooth hypersurface in R 3 and 
that Ug~(x) belongs to the uniform local Sobolev space H^Vq), while u^(x) belongs 
to the Sobolev space H s (T>q), for some fixed s > 10. Assume also that there is a 
function a(a) € H s (Sq) so that (2.8) -(2.9) hold, and the compatibility conditions 
up to order s — 1 are satisfied on So by the initial data, together with the entropy 
condition: 

^ • n(a) + c(p+, S+) < a(a) < • n(a) + cfo , 5 "), (2.11) 

and the Majda stability condition: 

<P s(p + ,S+) \ p(p+,S+) -p(p-,S~) 

P+6+ +l ) p+c*(p+,S+) < [ ' 

Then there exists T > and a C 2 -hypersurface S(t) together with C 1 -functions 
u ± (t,x) defined for t G [0,T] so i/iai u(£, x) = u ± (t, x), (i,x) e I? ± , is the dis- 
continuous shock-front solution of the Cauchy problem (1.4) and (2.10) satisfying 
(2.6)-(2.7). 

Condition in (2.12) is always satisfied for shocks of any strength for polytropic 
gas with 7 > 1 and for sufficiently weak shocks for general equation of state. In 
Theorem 2.5, c = c(p, S) is the sonic speed. The Sobolev space H^iV^) is defined 
as follows: A vector function u is in H*^ provided that there exists some r > 
such that 

x — y 

max ||w ry u||ij-s < oo with cj ry (x)=u;( ), 



where to e C^°(M d ) is a function so that w(x) > 0, w(x) = 1 when |x| < |, and 
w(x) = when |x| > 1. 

The compatibility conditions in Theorem 2.5 are defined in [156] and needed 
in order to avoid the formation of discontinuities in higher derivatives along other 
characteristic surfaces emanating from So- Once the main condition in (2.8) is 
satisfied, the compatibility conditions are automatically guaranteed for a wide class 
of initial functions. Further studies on the local existence and stability of shock- 
front solutions can be found in Majda [156, 157]. The existence of shock- front 
solutions whose lifespan is uniform with respect to the shock strength was obtained 
in Metivier [162]. Also see Blokhin-Trokhinin [20] for further discussions. 

The idea of the proof is similar to that for Theorem 2.2 by using the existence of 
a strictly convex entropy and the symmetrization of (1.1), but the technical details 
arc quite different due to the unusual features of the problem in Theorem 2.5. For 
more details, see [156]. 

The Navier-Stokes regularization of M-D shocks for the Euler equations has 
been established in Gues-Metiver- Williams- Zumbrun [120]. The local existence of 
rarefaction wave front-solutions for M-D hyperbolic systems of conservation laws 
has also been established in Alinhac [1]. Also see Benzoni-Gavage-Serre [14]. 

2.2. Hyperbolicity. There are many examples of m x m hyperbolic systems of 
conservation laws for d = 2 which are strictly hyperbolic; that is, they have simple 
characteristics. However, for d — 3, the situation is different. The following result 
for the case m = 2(mod4) is due to Lax [138] and for the more general case 
to = ±2,±3,±4(TOod8) due to Friedland-Robin-Sylvester [107]. 
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Theorem 2.6. Let A, B 7 and C be three mxm matrices such that aA + /3B + 7C 
has real eigenvalues for any real a, j3, and 7. Ifm = ±2, ±3, ±4 (mod 8), then there 
exist ao, A), arl ^ 7o witft a o + /^o +7o 7^ swc/i i/iai aoA + /3oB + 70C is degenerate, 
that is, there are at least two eigenvalues which coincide. 

This implies that, for d — 3, there are no strictly hyperbolic systems if to = 
±2, ±3, ±4 (mod 8). Such multiple characteristics influence the propagation of sin- 
gularities. 

Consider the isentropic Euler equations (1.9). When d = 2 and to — 3, the 
system is strictly hyperbolic with three real eigenvalues A_ < Ao < A + : 

A = n\V\ + n 2 v 2 , A± = mvi + n 2 v 2 ± c(p), p>0. 

Strict hyperbolicity fails at the vacuum states p — 0. However, when d = 3 and 
to = 4, the system is no longer strictly hyperbolic even when p > since the 
eigenvalue Ao = n\vi + n 2 v 2 + 713^3 has double multiplicity, although the other 
eigenvalues A± — n\v\ + n 2 v 2 + n^v^, ± c{p) are simple when p > 0. 

Consider the adiabatic Euler equations (1.4). When d = 2 and to = 4, the 
system is nonstrictly hyperbolic, since the eigenvalue Ao = + n 2 v 2 has double 
multiplicity, although A± = niVi+n 2 v 2 ±c(p, S) are simple when p > 0. When d = 3 
and to = 5, the system is again nonstrictly hyperbolic, since the eigenvalue A = 
niVi+n 2 v 2 +n^V3 has triple multiplicity, although A± = niWi+n 2 ^2+«3W3±c(p, 5) 
are simple when p > 0. 

2.3. Genuine Nonlinearity. The j t/l -characteristic field of system (1.1) in T> is 
called genuinely nonlinear if, for each fixed n e 5 fd_1 , the j t/l -cigcnvalue Aj(u;n) 
and the corresponding eigenvector r,-(u;n) determined by (Vf(u) • n)rj(u;n) = 
Xj (u; n)rj (u; n) satisfy 

V u Aj (u; n) • rj (u; n) 7^ for any u e D.neS". (2.13) 

The j t/l -charactcristic field of (1.1) in V is called linearly degenerate if 

V u Aj(u; n) • 1^(11; n) = for any u e V. (2.14) 

Then any scalar quasilinear conservation law in M. d ,d > 2, is never genuinely 
nonlinear in all directions. This is because, in this case, A(u;n) = f"(u) ■ n and 
r(u; n) = 1, and X'(u; n)r(u; n) = f"(u) • n, which is impossible to make this never 
equal to zero in all directions. A M-D version of genuine nonlinearity for scalar 
conservation laws is that the set 

{u : t + f'(u) • n = 0} has zero Lebesgue measure for any (r, n) G S , 

which is a generalization of (2.13). Under this generalized nonlinearity, the fol- 
lowing have been established: (i) Solution operators are compact as t > in 
Lions-Perthame-Tadmor [152] (also see [53, 193]); (ii) Decay of periodic solutions 
(Chen-Frid [52]; also see Engquist-E [102]); (iii) Strong initial and boundary traces 
of entropy solutions (Chen-Rasclc [60], Vasseur [202]; also see Panov [173]); (iv) 
SV-like structure of L°° entropy solutions (De Lellis-Otto-Westdickenberg [92]). 
Furthermore, we have 

Theorem 2.7 (Lax [138]). Every real, strictly hyperbolic quasilinear system for 
to = 2k, k > 1 odd, and d = 2 is linearly degenerate in some direction. 
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Quite often, linear degeneracy results from the loss of strict hyperbolicity. For 
example, even in the 1-D case, if there exists i ^ j such that A,(u) = Aj(u) for 
all u e K, then Boillat [22] proved that the i th - and ^-characteristic fields are 
linearly degenerate. 

For the isentropic Euler equations (1.9) with d — 2 and m = 3 for polytropic 
gases with the eigenvalues: 

A = nivi + n 2 v 2 , A± = n^vi + n 2 v 2 ± c(p), 

and the corresponding eigenvectors i"o and r±, we have VAo-ro = 0, which is linearly 
degenerate, and VA± • r± ^ for p e (0, oo), which are genuinely nonlinear. For 
the adiabatic Euler equations (1.4) with d = 2 and m = 4 for polytropic gases with 
the eigenvalues: 

A = nivi + n 2 v 2 , A± = n x v\ + n 2 v 2 ± c(p, S), 

and the corresponding eigenvectors r and r±, we have VAo-r = 0, which is linearly 
degenerate, and VA± • r± ^ for p, S e (0, oo), which are genuinely nonlinear. 

2.4. Singularities. For the 1-D case, singularities include the formation of shocks, 
contact discontinuities, and the development of vacuum states, at least for gas 
dynamics. For the M-D case, the situation is much more complicated: besides 
shocks, contact discontinuities, and vacuum states, singularities may include vor- 
ticity waves, focusing waves, concentration waves, complicated wave interactions, 
among others. Consider the Cauchy problem of the Euler equations in (1.4) in R 3 
for polytropic gases with smooth initial data: 

u| t =o = uo(x) with po(x) > 0, ( 2 -15) 

satisfying (po, mo, S*o)(x) = (p, 0, S) for |x| > L, where p > 0, S, and L arc 
constants. The equations in (1.4) possess a unique local C 1 -solution u(t,x) with 
p(t,x) > provided that the initial function (2.15) is sufficiently regular (Theorem 
2.2). The support of the smooth disturbance (po(x) — p, m (x), 5o(x) — S) propa- 
gates with speed at most c = \fp p {p 1 S) (the sound speed) that is, (p, m, 5*)(i, x) = 
(p, 0, S) if |x| > L + ct. Take p = p(p, S). Define 

P(t)= [ (p(t,x) 1/7 -p 1/7 )dx, F(t)= [ m(t,x)-xdx, 

which measure the entropy and the radial component of momentum. 

Theorem 2.8 (Sidcris [190]). Suppose that (p, m, S)(t,x) is a C 1 -solution of (1.4) 
and (2.15) for < t < T and 

P(0) > 0, F(0) > ^7rci 4 max( / 9 (x)). (2.16) 

Then the lifespan T of the C 1 -solution is finite. 

In particular, when po = P and So = S, P(0) = and the second condition in 
(2.16) holds if the initial velocity vo(x) = satisfies J\ X \ <L ^o( x ) • xdx > caL 4 . 

From this, one finds that the initial velocity must be supersonic in some region 
relative to c. The formation of a singularity (presumably a shock) is detected 
as the disturbance overtakes the wave-front forcing the front to propagate with 
supersonic speed. The formation of singularities occurs even without condition of 
largeness such as (2.16). For more details, see [190]. 
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In Christodoulou [79] , the relativistic Euler equations for a perfect fluid with an 
arbitrary equation of state have been analyzed. The initial function is imposed on 
a given spacelike hyperplane and is constant outside a compact set. Attention is 
restricted to the evolution of the solution within a region limited by two concentric 
spheres. Given a smooth solution, the geometry of the boundary of its domain of 
definition is studied, that is, the locus where shocks may form. Furthermore, under 
certain smallness assumptions on the size of the initial data, a remarkable and 
complete picture of the formation of shocks in R 3 has been obtained. In addition, 
sharp sufficient conditions on the initial data for the formation of shocks in the 
evolution have been identified, and sharp lower and upper bounds for the time of 
existence of a smooth solution have been derived. Also see Christodoulou [80] for 
the formation of black holes in general relativity. 

2.5. BV Bound. For 1-D strictly hyperbolic systems, Glimm's theorem [113] in- 
dicates that, as long as ||u ||w is sufficiently small, the solution u(t, x) satisfies 
the following stability estimate: 

||u(t,-)|| B v<C||uo||Bv (2-17) 

Even more strongly, for two solutions u(t, x) and w(i, x) obtained by either the 
Glimm scheme, wave-front tracking method, or vanishing viscosity method with 
small total variation, 

||u(t, •) - w(t, -)\\mm < c ||w(o, •) - v(o, -)IUi(R)- 

See Bianchini-Bressan [18] and Bressan [27]; also see Dafermos [91], Holden-Risebro 
[126], LeFloch [141], Liu- Yang [154], and the references cited therein. 

The recent great progress for entropy solutions to 1-D hyperbolic conservation 
laws based on BV estimates and trace theorems of BV fields naturally arises the 
expectation that a similar approach may also be effective for M-D hyperbolic sys- 
tems of conservation laws, that is, whether entropy solutions satisfy the relatively 
modest stability estimate (2.17). Unfortunately, this is not the case. Rauch [177] 
showed that the necessary condition for (2.17) to be held is 

Vf fe (u)Vf)(u) = Vfj(u)Vffc(u) for all k,l = 1,2,..., d. (2.18) 

The analysis above suggests that only systems in which the commutativity relation 
(2.18) holds offer some hope for treatment in the BV framework. This special 
case includes the scalar case m — 1 and the 1-D case d = 1. Beyond that, it 
contains very few systems of physical interest. An example is the system with 
fluxes: ffe(u) = (/>/j(|u| 2 )u, k = 1,2,..., d, which governs the flow of a fluid in 
an anisotropic porous medium. However, the recent study in Bressan [28] and 
Ambrosio-DeLellis [4] shows that, even in this case, the space BV is not a well- 
posed space (also cf. Jenssen [130]). Moreover, entropy solutions generally do not 
have even the relatively modest stability: ||u(f, •) — u||lp < C'pljuo — Uo||lp, p 7^ 2, 
based on the linear theory by Brenner [26] . 

In this regard, it is important to identify good analytical frameworks for the 
analysis of entropy solutions of M-D conservation laws (1.1), which are not in BV, 
or even in L p . A general framework is the space of divergence- measure fields, 
formulated recently in [54, 64, 66], which is based on entropy solutions satisfying 
the Lax entropy inequality (2.2). See §7 for more details. 

Another important notion of solutions is the notion of the measure- valued en- 
tropy solutions, based on the Young measure representation of a weak convergent 
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sequence (cf. [195, 12]), proposed by DiPerna [97]. An effort has been made to 
establish the existence of measure-valued solutions to the M-D Euler equations for 
compressible fluids; see Ganbo-Wcstdickenberg [111]. 

2.6. Nonuniqueness. Another important feature is the nonuniqueness of entropy 
solutions in general. In particular, De Lellis-Szekelyhidi [93] recently showed the 
following remarkable fact. 

Theorem 2.9. Let d > 2. Then, for any given function p = p(p) with p'{p) > 
when p > 0, there exist bounded initial functions (/>o,Vo) = (po, ^) with po( x ) > 
$o > f or which there exist infinitely many bounded solutions (p, m) of (1.9) with 
p> S for some 5 > 0, satisfying the energy identity in the sense of distributions: 

d t (pE)+V x - (m(E + p/p)) =0. (2.19) 

In fact, the same result also holds for the full Euler system, since the solutions 
constructed satisfy the energy equality so that there is no energy production at 
all. The main point for the result is that the solutions constructed contain only 
vortex sheets and vorticity waves which keep the energy identity (2.19) even for 
weak solutions in the sense of distributions, while the vortex sheets do not appear 
in the 1-D case. The construction of the infinitely many solutions is based on a 
variant of the Baire category method for differential inclusions. Therefore, for the 
uniqueness issue, we have to narrow down further the class of entropy solutions to 
single out physical relevant solutions, at least for the Euler equations. 

3. MULTIDIMENSIONAL SYSTEMS AND MODELS 

M-D problems are extremely rich and complicated. Some great developments 
have been made in the recent decades through strong and close interdisciplinary 
interactions and diverse approaches, including 

(i) Experimental data; 

(ii) Large and small scale computing by a search for effective numerical meth- 
ods; 

(iii) Asymptotic and qualitative modeling; 

(iv) Rigorous proofs for prototype problems and an understanding of the solu- 
tions. 

In some sense, the developments made by using approach (iv) are still behind those 
by using the other approaches (i)-(iii); however, most scientific problems are con- 
sidered to be solved satisfactorily only after approach (iv) is achieved. In this sec- 
tion, together with Sections 4-7, we present some samples of M-D systems/models 
and related important problems with emphasis on those prototypes that have been 
solved or may be expected to be solved rigorously at least for some cases. 

Since the M-D problems are so complicated in general, a natural strategy to 
attack these problems as a first step is to study (i) simpler nonlinear systems with 
strong physical motivations and (ii) special, concrete nonlinear physical problems. 
Meanwhile, extend the results and ideas from the first step to (i) study the Eu- 
ler equations in gas dynamics and elasticity; (ii) study more general problems; 
(iii) study nonlinear systems that the Euler equations are the main subsystem or 
describe the dynamics of macroscopic variables such as Navier- Stokes equations, 
MHD equations, combustion equations, Euler-Poisson equations, kinetic equations 
especially including the Boltzmann equation, among others. 
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Now we first focus on some samples of M-D systems and models. 

3.1. Unsteady Transonic Small Disturbance Equation. A simple model in 
transonic aerodynamics is the UTSD equation or so-called the 2-D inviscid Burgers 
equation (see Cole-Cook [81]): 

{ d t u + d x (ju 2 ) + d y v = 0, , . 

^ d y u — d x v = 0, ' ' 

or in the form of Zabolotskaya-Khokhlov equation [208] : 

d t (d t u + ud x u) + dyyU = 0. (3.2) 

The equations in (3.1) describe the potential flow field near the reflection point in 
weak shock reflection, which determines the leading order approximation of geomet- 
ric optical expansions; and it can also be used to formulate asymptotic equations 
for the transition from regular to Mach reflection for weak shocks. See Morawetz 
[165], Hunter [128], and the references cited therein. Equation (3.2) also arises in 
many different situations; see [128, 199, 208]. 

3.2. Pressure- Gradient Equations and Nonlinear Wave Equations. The 

inviscid fluid motions are driven mainly by the pressure gradient and the fluid 
convection (i.e., transport). As for modeling, it is natural to study hrst the effect 
of the two driving factors separately. 

Separating the pressure gradient from the Euler equations (1.4), we first have 
the pressure-gradient system: 

d t p = 0, d t (pv) + \7 xP = 0, a t (H+pV x -v = 0. (3.3) 

We may choose p = 1. Setting p = (7 — 1)P and t = and eliminating the 

velocity v, we obtain the following nonlinear wave equation for P: 

d ss (lnP)-A x P = 0. (3.4) 

Although system (3.3) is obtained from the splitting idea, it is a good approximation 
to the full Euler equations, especially when the velocity v is small and the adiabatic 
exponent 7 > 1 is large. See Zheng [215]. 

Another related model is the following nonlinear wave equation proposed by 
Canic-Keyfitz-Kim [32]: 

d tt p - A xP (p) = 0, (3.5) 
where p — p{p) is the pressure-density relation for fluids. Equation (3.5) is obtained 
from (1.9) by neglecting the inertial terms, i.e., the quadratic terms in m. This 
yields the following system: 

d t p + V x • m = 0, d t m + V x p(p) = 0, (3.6) 

which leads to (3.5) by eliminating m in the system. 

3.3. Pressureless Euler Equations. With the pressure-gradient equations (3.3), 
the convection (i.e., transport) part of fluid flow forms the pressureless Euler equa- 
tions: 

dtp + V x • Gov) = 0, 

9 t (pv) + V x • (pv ® v) = 0, (3.7) 
d t (pE)+V x - (pEw) =0. 
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This system also models the motion of free particles which stick under collision; see 
Brenier-Grenier [25], E-Rykov-Sinai [99], and Zeldovich [210]. In general, solutions 
of (3.7) become measure solutions. 

System (3.7) has been analyzed extensively; cf. [23, 25, 59, 99, 118, 127, 142, 
188, 204] and the references cited therein. In particular, the existence of measure 
solutions of the Riemann problem was first presented in Bouchut [23] for the 1-D 
case. It has been shown that 5-shocks and vacuum states do occur in the Riemann 
solutions even in the 1-D case. Since the two eigenvalues of the transport equa- 
tions coincide, the occurrence as t > can be regarded as a result of resonance 
between the two characteristic fields. Such phenomena can also be regarded as the 
phenomena of concentration and cavitation in solutions as the pressure vanishes. 
It has been rigorously shown in Chen-Liu [59] for 7 > 1 and Li [142] for 7 = 1 that, 
as the pressure vanishes, any two-shock Riemann solution to the Euler equations 
tends to a <5-shock solution to (3.7) and the intermediate densities between the two 
shocks tend to a weighted 5- measure that forms the <5-shock. By contrast, any two- 
rarefaction-wave Riemann solution of the Euler equations has been shown in [59] 
to tend to a two-contact-discontinuity solution to (3.7), whose intermediate state 
between the two contact discontinuities is a vacuum state, even when the initial 
density stays away from the vacuum. 

3.4. Incompressible Euler Equations. The incompressible Euler equations take 
the form: 

d t v + V • (v ® v) + Vp = 0, , . 

V-v = (3 - 8) 
for the density-independent case, and 

dtp + v • (pv) = 0, 

9 t (pv) + V-(pv®v) + Vp = 0, (3.9) 
V-v = 

for the density-dependent case, where p should be regarded as an unknown. These 
systems can be obtained by formal asymptotics for low Mach number expansions 
from the Euler equations (1.4). For more details, see Chorin [78], Constantin [82], 
Hoff [125], Lions [150], Lions-Masmoudi [151], Majda-Bertozzi [158], and the refer- 
ences cited therein. 

System (3.8) or (3.9) is an excellent model to capture M-D vorticity waves by 
ignoring the shocks in fluid flow, while system (1.11) is an excellent model to capture 
M-D shocks by ignoring the vorticity waves. 

3.5. Euler Equations in Nonlinear Elastodynamics. The equations of nonlin- 
ear elastodynamics provide another excellent example of the rich special structure 
one encounters when dealing with hyperbolic conservation laws. In M 3 , the state 
vector is (v, F), where v € K 3 is the velocity vector and F is the 3x3 matrix- valued 
deformation gradient constrained by the requirement that dctF > 0. The system 
of conservation laws, which express the integrability conditions between v and F 
and the balance of linear momentum, reads 

d t F ia - d Xa Vi = 0, i, a =1,2, 3, 

dtv j -T,l=id Xe S j0 (F)=O, j = 1,2,3. [6AU) 

The symbol S stands for the Piola-Kirchhoff stress tensor, which is determined by 
the (scalar-valued) strain energy function cr(F): Sjp(F) = 9 gp F ^ ■ System (3.10) is 
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hyperbolic if and only if, for any vectors (,n£ S 3 , 

d 2 a(F) 



E 



dF ia dF jB 

l<i,j,a,P<3 31 



£i£jn a np > 0. (3.11) 



System (3.10) is endowed with an entropy pair: 

7 ? = ( 7(F) + i|v| 2 , q a = - VjSjaiF). 

l<j<3 

However, the laws of physics do not allow <r(F), and thereby r/, to be convex func- 
tions. Indeed, the convexity of a would violate the principle of material frame 
indifference <x(OF) = u(F) for all O G 50(3) and would also be incompatible with 
the natural requirement that er(F) — > oo as det F 4- or det F f oo (see Dafermos 
[89]). 

The failure of a to be convex is also the main source of complication in elasto- 
statics, where one is seeking to determine equilibrium configurations of the body by 
minimizing the total strain energy J a(F)dx. The following alternative conditions, 
weaker than convexity and physically reasonable, are relevant in that context: 

(i) Poly convexity in the sense of Ball [11]: cr(F) = g(F, F*, det F), where 
F* = (detF)F -1 is the adjugate of F (the matrix of cofactors of F) and 
g(F, G, w) is a convex function of 19 variables; 

(ii) Quasiconvexity in the sense of Morrey [167]; 
(hi) Rank-one convexity, expressed by (3.11). 

It is known that convexity =>■ polyconvexity =>■ quasiconvexity rank-one convex- 
ity, however, none of the converse statements is generally valid. 

It is important to investigate the relevance of the above conditions in elastody- 
namics. A first start was made in Dafermos [89] where it was shown that rank-one 
convexity suffices for the local existence of classical solutions, quasiconvexity yields 
the uniqueness of classical solutions in the class of entropy solutions in L°°, and 
polyconvexity renders the system symmetrizable (also see [176]). To achieve this 
for polyconvexity, one of the main ideas is to enlarge system (3.10) into a large, 
albeit equivalent, system for the new state vector (v, F, F*,w) with w — detF: 

d t w= $*»(0<)> (3-12) 

l<Q,i<3 

dtF* k = ^2 dx*{eap*,eijkFji3Vi), 1,k = 1,2,3, (3.13) 

l<a,0,i,j<3 

where e Q( g 7 and e^fe denote the standard permutation symbols. Then the enlarged 
system (3.10) and (3.12)-(3.13) with 21 equations is endowed a uniformly convex 
entropy n = ct(F, F*, w) + ^|v| 2 so that the local existence of classical solutions and 
the stability of Lipschitz solutions may be inferred directly from Theorem 2.3. See 
Dafermos [91], Demoulini-Stuart-Tzavaras [94], and Qin [176] for more details. 

3.6. Born-Infeld System in Electromagnetism. The Born-Infeld system is a 
nonlinear version of the Maxwell equations: 

d t B + cnv\(^) =0, ^£>-curl(^)=0, (3.14) 
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where W : M 3 x K 3 — »■ K is the given energy density. The Born-Infeld model 
corresponds to the special case 

W B i(B,D) = ^T+W+W+W- 

When W is strongly convex (i.e., D 2 W > 0), system (3.15) is endowed with a 
strictly convex entropy. However, Wbi is not convex for a large enough field. As 
in §3.5, the Born- Infold model is enlarged from 6 to 10 equations in Brenier [24], 
by an adjunction of the conservation laws satisfied by P := B x D and W so that 
the augmented system turns out to be a set of conservation laws in the unknowns 
(h, B,J),F)eMxM 3 xl 3 xM 3 , which is in the physical region: 

{(h,B,D,P) : P = D x B,h= ^1 + \B\ 2 + \D\ 2 + \P\ 2 > 0}. 

Then the enlarged system is 

' dth + dWP = 0, 

d t B + curl( Px f+ D ) = 0, 
< $£> + cur l(£x£=*) =0 , ( 3 - 15 ) 
d t P + Div ( i^-BSB-Psg-r ) = o, 

which is endowed with a strongly convex entropy, where / is the 3x3 identity 
matrix. Also see Serre [185] for another enlarged system consisting of 9 scalar 
evolution equations in 9 unknowns (B,D,P), where P stands for the relaxation of 
the expression D x B. 



3.7. Lax Systems. Let /(u) be an analytic function of a single complex variable 
u = u + vi. We impose on the complex- valued function u = u(t, z), z = x + yi, and 
the real variable t the following nonlinear PDE: 

d t u + d z f(u) = 0, (3.16) 

where the bar denotes the complex conjugate and d z — \{d x — id y ). We may 
express this equation in terms of the real and imaginary parts of u and |/(u) = 
a(u, v) + b(u,v)i. Then (3.16) gives 

d t u + d x a(u,v) + d y b(u,v) = 0, 

dtv — d x b(u, v) + d y a(u, v) = 0. ^ ' ' 

In particular, when /(u) = u 2 = u 2 + v 2 + 2uvi, system (3.16) is called the complex 
Burger equation, which becomes 

( d t u+±d x (u 2 +v 2 )+d y (uv)=0, 

\ d t v - d x (uv) + \d y {u 2 + v 2 ) = o. ^- 18 J 

System (3.17) is a symmetric hyperbolic system of conservation laws with a strictly 
convex entropy r](u, v) = u 2 + v 2 ; see Lax [139] for more details. For the 1-D case, 
this system is an archetype of hyperbolic conservation laws with umbilic degener- 
acy, which has been analyzed in Chen-Kan [56], Schacffer-Shearcr [181], and the 
references cited therein. 
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3.8. Gauss-Codazzi System for Isometric Embedding. A fundamental prob- 
lem in differential geometry is to characterize intrinsic metrics on a 2-D Ricmannian 
manifold M 2 which can be realized as isometric immersions into R 3 (cf. Yau [207]; 
also see [124]). For this, it suffices to solve the Gauss-Codazzi system, which can 
be written as (cf. [62, 124]) 

( d x M -d v L = T^L-2T { ^M + T^N, 



d x N - d v M = -Y { })L + 2lil ) M - tWn, 



(3.19) 

-22 ^ "•" * L 12 1V1 L 11 A 

and 

LN-M 2 = K. (3.20) 

Here K(x, y) is the given Gauss curvature, and is the given Christoffel symbol. 
We now follow Chen-Slemrod-Wang [62] to present the fluid dynamical formulation 
of system (3.19)-(3.20). Set L = pv 2 + p,M = -puv,N = pu 2 + p, and set 
q 2 = u 2 + v 2 as usual. Then the Codazzi equations (3.19) become the familiar 
balance laws of momentum: 

f d x (puv) + d y (pv 2 +p) = -T<§(pv 2 +p)- 2Y { £puv - T^(pu 2 + p), ^ 

\ d x (pu 2 +p) + d y (puv) = -T$(pv 2 +p) - 2T[^puv - T { ^{pu 2 +p), 

and the Gauss equation (3.20) becomes ppq 2 + p 2 = K. We choose p to be the 
Chaplygin gas- type to allow the Gauss curvature K to change sign: p = —K Then 
we obtain the "Bernoulli" relation: 

or p = -y / q 2 + K. (3.22) 



In general, system (3.21)-(3.22) for unknown (u, v) is of mixed hyperbolic-elliptic 
type determined by the sign of Gauss curvature (hyperbolic when K < and elliptic 
when K > 0). 

For the Gauss-Codazzi- Ricci equations for isometric embedding of higher dimen- 
sional Riemannian manifolds, see Chen-Slemrod-Wang [63]. 



4. Multidimensional Steady Supersonic Problems 

M-D steady problems for the Euler equations are fundamental in fluid mechan- 
ics. In particular, understanding of these problems helps us to understand the 
asymptotic behavior of evolution solutions for large time, especially global attrac- 
tors. One of the excellent sources of steady problems is Courant-Friedrichs's book 
[84] . In this section we first discuss some of recent developments in the analysis of 
M-D steady supersonic problems. The M-D steady Euler flows are governed by 

V x • (pv) = 0, 

V x (pv®v) + V x p = 0, (4.1) 
V x -(p(£+2)v)=0, 

where v is the velocity, E is the total energy, and the constitutive relations among 
the thermodynamical variables p,p 7 e 7 8 7 and S are determined by (1.5)— (1.7) . 

For the barotropic (isentropic or isothermal) case, p = p(p) is determined by 
(1.10) with 7 > 1, and then the first d+1 equations in (4.1) form a self-contained 
system, the Euler system for steady barotropic fluids. 
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System (4.1) governing a supersonic flow (i.e., |v| 2 > c 2 ) has all real eigenvalues 
and is hyperbolic, while system (4.1) governing a subsonic flow (i.e., |v| 2 < c 2 ) has 
complex eigenvalues and is elliptic-hyperbolic mixed and composite. 

4.1. Wedge Problems Involving Supersonic Shock- Fronts. The analysis of 
2-D steady supersonic flows past wedges whose vertex angles are less than the criti- 
cal angle can date back to the 1940s since the stability of such flows is fundamental 
in applications (cf. [84, 205]). Local solutions around the wedge vertex were first 
constructed in Gu [119], Li [146], and Schaeffer [180]. Also see Zhang [212] and the 
references cited therein for global potential solutions when the wedges are a small 
perturbation of the straight-sided wedge. For the wedge problem, when the vertex 
angle is suitably large, the flow contains a large shock-front and, for this case, the 
full Euler equations (4.1) are required to describe the physical flow. When a wedge 
is straight and its vertex angle is less than the critical angle ui c , there exists a su- 
personic shock-front emanating from the wedge vertex so that the constant states 
on both sides of the shock are supersonic; the critical angle condition is necessary 
and sufficient for the existence of the supersonic shock (also see [69, 84]). 

Consider 2-D steady supersonic Euler flows past a 2-D Lipschitz curved wedge 
\ x 2\ < 5(^1)^1 > 0, with g e Lip(R + ) and g(0) = 0, whose vertex angle uo '■= 
arctan(g'(0+)) is less than the critical angle u) c , along which TV{g'(-); M+} < e for 
some constant e > 0. Denote 

fl := {x : x 2 > g( Xl ), Xl > 0}, r := {x : x 2 = g( Xl ), x x > 0}, 

and n(xi±) = ( a (xi±)s) ^ outer normal vectors to T at points ii±, 

respectively. The uniform upstream flow u = (/?_, 0, p-EJ) satisfies u_ > 

c(/9_, S-) so that a strong supersonic shock- front emanates from the wedge vertex. 
Since the problem is symmetric with respect to the X2-axis, the wedge problem can 
be formulated into the following problem of initial-boundary value type for system 
(4.1) in O: 

Cauchy Condition: u\ Xl =v = u_; (4.2) 

Boundary Condition : v • n = on T. (4.3) 

In Chen-Zhang-Zhu [69], it has been established that there exist £0 > and 
C > such that, when s < £0, there exists a pair of functions 

u = (p, v, pE) e BV(R; M 2 x R+ x E+), a E BV(R+:R) 

with x = Jq 1 o-(s)ds e Lip(M + ;M + ) such that u is a global entropy solution of 
problem (4.1)-(4.3) in tt with 

TV{u(xi, •) : [g{ Xl ), -00)} < CTV{g'(-)} for every Xl G R+, (4.4) 

and the strong shock- front x 2 = x{ x i) emanating from the wedge vertex is nonlin- 
early stable in structure. Furthermore, the global ^-stability of entropy solutions 
with respect to the incoming flow at X \ = in L 1 has also been established in 
Chen-Li [58]. This asserts that any supersonic shock- front for the wedge problem 
is nonlinearly Z^-stable with respect to the BV perturbation of the incoming flow 
and the wedge boundary. 

In order to achieve this, we have first developed an adaptation of the Glimm 
scheme whose mesh grids are designed to follow the slope of the Lipschitz wedge 
boundary so that the lateral Riemann building blocks contain only one shock or 
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rarefaction wave emanating from the mesh points on the boundary. Such a design 
makes the BV estimates more convenient for the Glimm approximate solutions. 
Then careful interaction estimates have been made. One of the essential estimates 
is the estimate of the strength 5\ of the reflected 1-waves in the interaction between 
the 4-strong shock-front and weak waves (a\, (3 2 , 03, fii), that is, 

tfi=ai+Jr a 04 + O(l)KI(|&| + |&|) with \K sl \<l. 

The second essential estimate is the interaction estimate between the wedge bound- 
ary and weak waves. Based on the construction of the approximation solutions and 
interaction estimates, we have successfully identified a Glimm- type functional to 
incorporate the curved wedge boundary and the strong shock-front naturally and 
to trace the interactions not only between the wedge boundary and weak waves, but 
also between the strong shock-front and weak waves. With the aid of the important 
fact that l-Kgil < 1, we have showed that the identified Glimm functional mono- 
tonically decreases in the flow direction. Another essential estimate is to trace the 
approximate strong shock-fronts in order to establish the nonlinear stability and 
asymptotic behavior of the strong shock-front emanating from the wedge vertex 
under the BV wedge perturbation. 

For the 3-D cone problem, the nonlinear stability of a self-similar 3-D full gas 
flow past an infinite cone is another important problem. See Lien-Liu [140] for the 
cones with small vertex angle. Also see [71, 77] for the construction of piecewise 
smooth potential flows under smooth perturbation of the straight-sided cone. 

In Elling-Liu [101], an evidence has been provided that the steady supersonic 
weak shock solution is dynamically stable, in the sense that it describes the long- 
time behavior of an unsteady flow. 

4.2. Stability of Supersonic Vortex Sheets. Another natural problem is the 
stability of supersonic vortex sheets above the Lipschitz wall x 2 — g(xi),xi > 0, 
with 

g e Lip{R + ;R), g(0) = g'{0+) = 0, lim arctan( 5 '(xi+)) = 0, 

and g' <G £?V(R+;R) such that TV{g'(-)} < e for some e > 0. Denote again 
Q = {x : x 2 > g{xi), x\ > 0} and F = {x : x 2 = g{x\), x\ > 0}. The upstream 
flow consists of one supersonic straight vortex sheet x 2 = yo > and two constant 
vectors u = (p , p u ,0, p E ) when x 2 > yo > and Ui = (pi, pim, 0, pi-Ei) 
when < x 2 < yo satisfying u\ > uq > and Ui > c(pi, Si) for i = 0, 1. Then 
the vortex sheet problem can be formulated into the following problem of initial- 
boundary value type for system (4.1): 

Cauchy Condition: uL 1=0 = / U °' < X2 < y °' (4.5) 

y 1X1 \ ui, x 2 >y ; y > 

Slip Boundary Condition : v • n = on T. (4.6) 

It has been proved that steady supersonic vortex sheets, as time-asymptotics, are 
stable in structure globally, even under the BV perturbation of the Lipschitz walls 
in Chcn-Zhang-Zhu [70]. The result indicates that the strong supersonic vortex 
sheets are nonlinearly stable in structure globally under the BV perturbation of 
the Lipschitz wall, although there may be weak shocks and supersonic vortex sheets 
away from the strong vortex sheet. In order to establish this theorem, as in §4.1, 
we first developed an adaption of the Glimm scheme whose mesh grids are designed 
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to follow the slope of the Lipschitz boundary. For this case, one of the essential 
estimates is the estimate of the strength 5i of the reflected 1-wave in the interaction 
between the 4-weak wave a± and the strong vortex sheet from below is less than 
one: 8\ = K m a^ and l-f^oil < 1- The second new essential estimate is the estimate 
of the strength £4 of the reflected 4-wave in the interaction between the 1-weak 
wave P\ and the strong vortex sheet from above is also less than one: <5 4 = Ku(3i 
and \Kn\ < 1. Another essential estimate is to trace the approximate supersonic 
vortex sheets under the BV boundary perturbation. For more details, see Chen- 
Zhang-Zhu [70]. The nonlinear ^-stability of entropy solutions with respect to the 
incoming flow at x\ = is under current investigation. 

5. Multidimensional Steady Transonic Problems 

In this section we discuss another important class of M-D steady problems: tran- 
sonic problems. In the last decade, a program has been initiated on the existence 
and stability of M-D transonic shock-fronts, and some new analytical approaches in- 
cluding techniques, methods, and ideas have been developed. For clear presentation, 
we focus mainly on the celebrated steady potential flow equation of aerodynamics 
for the velocity potential cp : O C R d — > M, which is a second-order nonlinear PDE 
of mixed elliptic-hyperbolic type: 

V x • (p(|V x </>| 2 )V x </>) = 0, xeOcM d , (5.1) 

where the density p(q 2 ) by scaling is 

ptf)=(l-l^±f)& ( 5 . 2 ) 

with adiabatic exponent 7 > 1. Equation (5.1) is elliptic at V x </? with |V X (^| = q if 
p(q 2 ) + 2q 2 p'{q 2 ) > 0, which is equivalent to 

q < q, := ^2/(7+1), 

i.e., the flow is subsonic. Equation (5.1) is hyperbolic if p(q 2 ) + 2q 2 p'(q 2 ) < 0, i.e., 
q > (?*, that is, the flow is supersonic. 

Let fl + and 0~ be open subsets of f2 such that fl + n f2~ = 0, f2+ U f2~ = Q, 
and S = dn+ n 0. Let ip g C^ij) n C a (IW) be a weak solution of (5.1), which 
satisfies |V x v| < Qcav '■= \/2/(7 — 1)j so that ^x<P experiences a jump across S 
that is a (d — 1)-D smooth surface. Set ^ — <£>|n±. Then ip satisfies the following 
Rankinc-Hugoniot conditions on S: 

ip+ = ip~, (5.3) 
p(|V x ^+| 2 )V x ^+ • n = p(|V x ^-| 2 )V x ^- • n. (5.4) 

Suppose that tp e C 1 (fi^) is a weak solution satisfying 

|V X (^| < in \V x <p\ > in fT, V x ^ • n| 5 > 0. (5.5) 

Then <p is a transonic shock-front solution with transonic shock-front S dividing 
n into the subsonic region f2 + and the supersonic region and satisfying the 
physical entropy condition (see Courant-Friedrichs [84]): 

p(|V x (^-| 2 ) <p(|V x ^+| 2 ) along S. (5.6) 

Note that equation (5.1) is elliptic in the subsonic region and hyperbolic in the 
supersonic region. 
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As an example, let (xi,x') be the coordinates in R d , where x\ G R and x' = 
(x 2 ,...,x d ) G Fix V G R d , and let <^ (x) := V -x, x G M d . If |V | G (0,g*) 

(resp. | Vo| G (<?*, (feat;))) then <po( x ) is a subsonic (resp. supersonic) solution in R d , 
and V = V x (^o is its velocity. 

Let Qq > and V G R d_1 be such that the vector Vq := (qq ,~V' ) satisfies 
|Vq~| > g*. Then there exists a unique > such that 

(1 - ^(koT + |V | 2 ))^V = (1 - ^(ko I 2 + |V | a ))^flb ■ (5-7) 

The entropy condition (5.6) implies < q$ . By denoting Vg" := (qo ,V ) an d 
defining the functions ip^(x) := ■ x on R d , then tp£ (resp. ip^) is a subsonic 
(resp. supersonic) solution. Furthermore, from (5.4) and (5.7), the function 

<A)(x) := min(^ -(x),^+(x)) = • x, x G fij := {±ari > 0}, (5.8) 

is a plane transonic shock- front solution in R d , Oq and Oq are respectively its 
supersonic and subsonic regions, and S = {x\ = 0} is a transonic shock-front. 
Note that, if V = 0, the velocities are orthogonal to the shock-front S and, if 
V 7^ 0, the velocities are not orthogonal to S. 

5.1. Transonic Shock-Front Problems in R d . Consider M-D perturbations of 
the uniform transonic shock-front solution (5.8) in R d with d > 3. Since it suffices to 
specify the supersonic perturbation ip~ only in a neighborhood of the unperturbed 
shock-front {x\ — 0}, we introduce domains f2 := (— l,oo) x R d_1 ,£2i := (—1, 1) x 
R d_1 . Note that we expect the subsonic region 0+ to be close to the half-space 

n+ = {x! > o}. 

Problem 5.1. Given a supersonic solution <p~(x.) of (5.1) in Qi, find a transonic 
shock-front solution </?(x) in such that 

f2 _ C fii, V 5 ( x ) = V _ ( x ) i n 

where ft~ := f2 \ il + and 51+ := {x G : |V x <^(x)| < q*}, and 

ip = ip~, d Xl <p = d Xl tp~ on {ari = -1}, (5.9) 

lim ||^-^||ci(n+\B s (o)) =0. (5.10) 

Condition (5.9) determines that the solution has supersonic upstream, while 
condition (5.10) especially determines that the uniform velocity state at infinity in 
the downstream direction is equal to the unperturbed downstream velocity state. 
The additional requirement in (5.10) that ip — > at infinity within fl + fixes 
the position of shock-front at infinity. This allows us to determine the solution of 
Problem 5.1 uniquely. In Chen-Feldman [47], we have employed the free boundary 
approach first developed for the potential flow equation (5.1) to solve the stability 
problem, Problem 5.1. This existence result can be extended to the case that the 
regularity of the steady perturbation ip~ is only C 1 ' 1 . It would be interesting to 
establish similar results for the full Euler equations (1.4). 

5.2. Nozzle Problems Involving Transonic Shock-Fronts. We now consider 
M-D transonic shock-fronts in the following infinite nozzle f2 with arbitrary smooth 
cross-sections: il = ^(A xK)(l {x\ > —1}, where A C R^" 1 is an open bounded 
connected set with a smooth boundary, and ^ : R d — > M. d is a smooth map, which 
is close to the identity map. For simplicity, we assume that d A is in 

C[*]+3,a and 
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I* — J||[d] + 3 a Rd < v for some a G (0, 1) and small a > 0, where [s] is the integer 
part of s, I : WL d -> M d is the identity map, and d t Q, := *(R x <9A) n {xi > -1}. 
For concreteness, we also assume that there exists L > 1 such that ^(x) = x for 
any x = (xi,x') with x\ > L. 

In the 2-D case, il = {(x\,X2) ■ Xi > — 1, 6~(x 2 ) < x 2 < b + (x 2 )}, where 
||fe ± — &^||4, q .r < <J and ^ = on [L, oo) for some constants 6^ satisfying 
&+ > fo^o- For the M-D case, the geometry of the nozzles is much richer. Note that 
our setup implies that <9f2 = d Q U d/fi with 

d t fl := *[(-oo,oo) x dA] n {(xi,x') : x x > -1}, 
d o := *((— oo, oo) x A) n {(ari, x') : x x = -1}. 
Then our transonic nozzle problem can be formulated as 

Problem 5.2: Transonic Nozzle Problem. Given the supersonic upstream flow 
at the entrance d $l: 

ip = <p~, ip Xl =ip~ on <9 fi, (5-11) 
the slip boundary condition on the nozzle boundary <9;fi: 

V x ip ■ n = on dift, (5.12) 

and the uniform subsonic flow condition at the infinite exit x\ = oo: 

\\<p(-) ~ qooXiWc^nnix^R}) -> as R -> oo for some q x G (0,q*), (5.13) 

find a M-D transonic flow of problem (5.1) and (5.11)— (5.13) in O. 

The standard local existence theory of smooth solutions for the initial-boundary 
value problem (5.11)— (5.12) for second-order quasilinear hyperbolic equations im- 
plies that, as a is sufficiently small, there exists a supersonic solution ip~ of (5.1) 
in i7 2 := { — 1 < X\ < 1}, which is a C l+1 perturbation of <£q = q^xf. For any 
ae(0,l], 

\W~ -Vo\\i,cc,n 2 <C o, 1 = 1,2, (5.14) 

for some constant Co > 0, and satisfies V(^>~ -n = on <9;f2 2 , provided that (ipj 
on <9 fi satisfying 

\\<P7 - %Xi\\h>+i + \\ip7 - %\\h°+'-i < ' = 1, 2, (5.15) 

for some integer s > d/2 + 1 and the compatibility conditions up to order s + 1, 
where the norm || • ||#s is the Sobolcv norm with H s = W s ' 2 . 

In Chcn-Fcldman [48], Problem 5.2 has been solved. More precisely, let G 
(q*,q C av) and q^ G (0,q*) satisfy (5.7), and let ipo be the transonic shock-front 
solution (5.8) with V = 0. Then there exist <ro > and C, depending only on d, 
a, 7, <7o~, A, and L, such that, for every a G (0, ao), any map H> introduced above, 
and any supersonic upstream flow (<fie ^e) 011 satisfying (5.15) with I = 1, 
there exists a solution (p G C 0,1 (f2) n C 1,a (f2+) satisfying 

n+fo>) = {xi > /(x')}, n-( v ) = {xi < j(x')}, 

Ik - q ± xi\\i, a ,si± < Co-, \\tp- qooXiWc^nnix^R}) -> as i? ->• oo, 
where <7oo G (0,q*) is the unique solution of the equation 

p{{q + f)q + = Q + ■= t^tt I P(|V X ^-| 2 + (C) 2 )C dW" -1 - (5-16) 
l A l Ja n 
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In addition, if the supersonic uniform flow (ip~,tp~) on d Vl satisfies (5.15) with 
/ = 2, then ip e C 2 ' a (fl + ) with \\ip — q + X\ ||2, a ,n + < Ca, and the solution with 
a transonic shock-front is unique and stable with respect to the nozzle boundary 
and the smooth supersonic upstream flow at the entrance. The techniques have 
been extended to solving the nozzle problem for the 2-D full Euler equations first 
in Chen-Chen- Feldman [43]. 

In the previous setting, the location of the transonic shock-front in the solution 
is not unique in general if we prescribe only the pressure at the nozzle exit, since 
the flat shock-front between uniform states can be translated along the flat nozzle 
which does not change the flow parameters at the entrance and the pressure at the 
exit. The analysis of the relation among the unique shock-front location, the flow 
parameters, and the geometry of the diverging nozzle has been made for various 
Euler systems in different dimensions; see Chen [73] and Lin- Yuan [153] for the 
2-D full Euler equations, Li-Xin-Yin [145] for the 2-D and 3-D axisymmetric Euler 
systems with narrow divergent nozzles, and Bae-Feldman [9] for perturbed diverging 
cone-shaped nozzle of arbitrary cross-section for the non-isentropic potential flow 
system in any dimension. 

A major open problem is the physical de Laval nozzle problem: Consider a nozzle 
which is flat (i.e., \I> = I) between —2 < x\ < — 1 and 1 < x\ < oo, and has some 
special geometry between — 1 < x\ < 1 to make the part { — 1 < x\ < 0} convergent 
and the part {0 < x\ < 1} divergent. Given certain incoming subsonic flow at 
x\ = —2, find a transonic flow containing transonic shock-fronts in the nozzle such 
that the downstream at x\ = oo is subsonic. A related reference is Kuz'min's book 
[135]. 

5.3. Wedge/Cone Problems Involving Transonic Shock-Fronts. The exis- 
tence and stability of transonic flows past infinite wedges or cones are further long- 
standing open transonic problems. Some progress has been made for the wedge case 
in 2-D steady flow in Chen- Fang [75], Fang [104], and Chen-Chen-Feldman [44]. In 
[104], it was proved that the transonic shock- front is conditionally stable under per- 
turbation of the upstream flow and/or the wedge boundary in some weak Sobolev 
norms. In [44], the existence and stability of transonic flows past the curved wedge 
have been established in the strong Holder norms for the full Euler equations. 

Conical flow (i.e. cylindrically symmetric flow with respect to an axis, say, the 
Zi-axis) in M 3 occurs in many physical situations (cf. [84]). Unlike the 2-D case, the 
governing equations for the 3-D conical case have a singularity at the cone vertex 
and the flow past the straight-sided cone is self-similar, but is no longer piecewise 
constant. These have resulted in additional difficulties for the stability problem. In 
Chen- Fang [46], we have developed techniques to handle the singular terms in the 
equations and the singularity of the solutions. Our main results indicate that the 
self-similar transonic shock-front is conditionally stable with respect to the conical 
perturbation of the cone boundary and the upstream flow in appropriate function 
spaces. That is, the transonic shock-front and downstream flow in our solutions arc 
close to the unperturbed self-similar transonic shock-front and downstream flow 
under the conical perturbation, and the slope of the shock-front asymptotically 
tends to the slope of the unperturbed self-similar shock-front at infinity. 

In order to achieve these results, we have first formulated the stability problem 
as a free boundary problem and have then introduced a coordinate transformation 
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to reduce the free boundary problem into a fixed boundary value problem for a 
singular nonlinear elliptic system. We have developed an iteration scheme that 
consists of two iteration mappings: one is for an iteration of approximate transonic 
shock-fronts; and the other is for an iteration of the corresponding boundary value 
problems for the singular nonlinear systems for given approximate shock-fronts. To 
ensure the well-definedness and contraction property of the iteration mappings, it is 
essential to establish the well-posedness for a corresponding singular linearized ellip- 
tic equation, especially the stability with respect to the coefficients of the equation, 
and to obtain the estimates of its solutions reflecting their singularity at the cone 
vertex and decay at infinity. The approach is to employ key features of the equa- 
tion, introduce appropriate solution spaces, and apply a Fredholm-type theorem in 
Maz'ya-Plamenevskivi [161] to establish the existence of solutions by showing the 
uniqueness in the solution spaces. Also see Cui-Yin [87] for related results. 

5.4. Airfoil/Obstacle Problems: Subsonic Flows past an Airfoil or an Ob- 
stacle. The Euler equations for potential flows (5.1)-(5.2) in M 2 can be rewritten 
as 



with (u, v) = V x <p and q = \V x <p\, where p is again determined by the Bernoulli re- 
lation (5.1). An important problem is subsonic flows past an airfoil or an obstacle. 
Shiffman [189], Bers [15], Finn-Gilbarg [106], and Dong [98] studied subsonic (ellip- 
tic) solutions of (5.1) outside an obstacle when the upstream flows are sufficiently 
subsonic; also see Chen-Dafermos-Slemrod-Wang [45] via compensated compactness 
argument. Morawetz in [166] first showed that the flows of (5.17) past an obstacle 
may contain transonic shocks in general. A further problem is to construct global 
entropy solutions of the airfoil problem (see [163, 110]). 

In Chen-Slemrod-Wang [61] , we have introduced the usual flow angle 9 = tan -1 (^) 
and written the irrotationality and mass conservation equation as an artificially vis- 
cous problem: 



where a(p) is suitably chosen, and appropriate boundary conditions are imposed 
for this regularized "viscous" problem. The crucial new discovery is that a uni- 
formly L°° bound in q 6 can be obtained when 1 < 7 < 3 which uniformly prevents 
cavitation. However, in this formulation, a uniform bound in the flow angle 9 £ and 
a uniform lower bound in q 6 in any fixed region disjoint from the profile must be as- 
sumed apriori. By making further careful energy estimates, Morawetz's argument 
[163] then applies, and the strong convergence in Ll oc (Q) of our approximating 
sequence is achieved. 

An important open problem is whether two remaining conditions on (q £ 7 9 £ ) can 
be removed. See Chen-Slemrod-Wang [61] for more details. 




(5.17) 




d Xl (p s u s )+d X2 (p s v s ) = eV • (<j( P e )Vp £ ) 



(5.18) 



5.5. Nonlinear Approaches. We now discuss several nonlinear approaches to 
deal with steady transonic problems. 
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5.5.1. Free Boundary Approaches. We first describe two of the free boundary ap- 
proaches for Problems 5.1-5.2, developed in [47, 48]. 

Free Boundary Problems. The transonic shock-front problems can be formu- 
lated into a one-phase free boundary problem for a nonlinear elliptic PDE: Given 
<p~ e C 1,a (0), find a function ip that is continuous in £1 and satisfies 



equation (5.1), the cllipticity condition in the non-coincidence set fi + = {<p < <p~}, 
the free boundary condition (5.4) on the boundary S = dfl + n f2, as well as the 
prescribed conditions on the fixed boundary <9f2 and at infinity. These conditions 
are different in different problems, for example, conditions (5.9)-(5.10) for Problem 
5.1 and (5.11) (5.13) for Problem 5.2. 

The free boundary is the location of the shock-front, and the free boundary con- 
ditions (5.3)-(5.4) are the Rankine-Hugoniot conditions. Note that condition (5.19) 
is motivated by the similar property (5.8) of unperturbed shock-fronts; and (5.19), 
locally on the shock- front, is equivalent to the entropy condition (5.6). Condition 
(5.19) transforms the transonic shock-front problem, in which the subsonic region 
£! + is determined by the gradient condition |V x <p(x)| < q*, into a free bound- 
ary problem in which tt + is the non-coincidence set. In order to solve this free 
boundary problem, equation (5.1) is modified to be uniformly elliptic and then the 
free boundary condition (5.4) is correspondingly modified. The problem becomes 
a one-phase free boundary problem for the uniformly elliptic equation which we 
can solve. Since ip~ is a small ^'"-perturbation of tp^ , the solution ip of the free 
boundary problem is shown to be a small ^'"-perturbation of the given subsonic 
shock- front solution tp^ in + . In particular, the gradient estimate implies that <p 
in fact satisfies the original free boundary problem, hence the transonic shock-front 
problem, Problem 5.1 (Problem 5.2, respectively). 

The modified free boundary problem does not directly fit into the variational 
framework of Alt-Caffarelli [2] and Alt-Caffarclli- Friedman [3], and the regulariza- 
tion framework of Berestycki-Caffarelli-Nirenberg [17]. Also, the nonlinearity of the 
free boundary problem makes it difficult to apply the Harnack inequality approach 
of Caffarelli [30]. In particular, a boundary comparison principle for positive solu- 
tions of nonlinear elliptic equations in Lipschitz domains is not available yet for the 
equations that are not homogeneous with respect to (V x u, V x u, u), which is our 
case. 

Iteration Approach. The first approach we developed in Chen-Fcldman [47, 
48] for the steady potential flow equation is an iteration scheme based on the non- 
degeneracy of the free boundary condition: the jump of the normal derivative of a 
solution across the free boundary has a strictly positive lower bound. Our iteration 
process is as follows. Suppose that the domain is given so that Sk ■= dQ^ \ <9£1 
is C l a . Consider the oblique derivative problem in obtained by rewriting the 
(modified) equation (5.1) and free boundary condition (5.4) in terms of the function 
u := <p — (£>q~. Then the problem has the following form: 



tp<ip 



in f2, 



(5.19) 



div x A(x, V x u) = F(x) 
A(x, V x u) • n = G(x, n) 



in fl+ := {u > 0}, 
on S := \ dfl, 



(5.20) 



plus the fixed boundary conditions on dQ^ fl dQ. and the conditions at infinity. 
The equation is quasilinear, uniformly elliptic, A(x, 0) = 0, while G(x, n) has a 



MULTIDIMENSIONAL CONSERVATION LAWS 



25 



certain structure. Let Uk € C 1,a (fijJ") be the solution of (5.20). Then \\uk\\ l a n + 
is estimated to be small if the perturbation is small, where appropriate weighted 
Holder norms are actually needed in the unbounded domains. The function ip^ := 
(p^+Uk from fijj" is extended to f2 so that the C la -norm of ipk — p^ in fi is controlled 
by \\uk\\ l a q+. For the next step, define := {x e ft : ^fc(x) < ip~(x.)}. Note 

that, since \\ipk — <Po \\i,a,n and \\<p~-tfo ||i,a,fi are small, we have |V x v?~|-|V x </5 fe | > 
S > in O, and this nondegeneracy implies that Sk+i ■= dil^ +1 \ dfi is C l a and 
its norm is estimated in terms of the data of the problem. 

The fixed point il + of this process determines a solution of the free boundary 
problem since the corresponding solution <p satisfies il + — {<p < p~} and the 
Rankine-Hugoniot condition (5.4) holds on S := dil + n f2. On the other hand, 
the elliptic estimates alone are not sufficient to get the existence of a fixed point, 
because the right-hand side of the boundary condition in problem (5.20) depends on 
the unit normal n of the free boundary. One way is to require the orthogonality of 
the flat shock-fronts so that p(\V x (p^ \ 2 )\/ x p^ = p(\V x <Pq | 2 )V x (/?q in fi to obtain 
better estimates for the iteration and prove the existence of a fixed point. For more 
details, see Chen-Feldman [47, 48]. 

Partial Hodograph Approach. The second approach we have developed in 
[47] is a partial hodograph procedure, with which we can handle the existence and 
stability of M-D transonic shock-fronts that are not nearly orthogonal to the flow 
direction. One of the main ingredients in this new approach is to employ a partial 
hodograph transform to reduce the free boundary problem to a conormal bound- 
ary value problem for the corresponding nonlinear second-order elliptic equation 
of divergence form and then to develop techniques to solve the conormal bound- 
ary value problem. To achieve this, the strategy is to construct first solutions in 
the intersection domains between the physical unbounded domain under consid- 
eration and a series of half balls with radius R, then make uniform estimates in 
R, and finally send R — > oo. It requires delicate apriori estimates to achieve this. 
A uniform bound in a weighted L°°-norm can be achieved by both employing a 
comparison principle and identifying a global function with the same decay rate as 
the fundamental solution of the elliptic equation with constant coefficients which 
controls the solutions. Then, by scaling arguments, the uniform estimates can be 
obtained in a weighted Holder norm for the solutions, which lead to the existence 
of a solution in the unbounded domain with some decay rate at infinity. For such 
decaying solutions, a comparison principle holds, which implies the uniqueness for 
the conormal problem. Finally, by the gradient estimate, the limit function can be 
shown to be a solution of the M-D transonic shock problem, and then the existence 
result can be extended to the case that the regularity of the steady perturbation 
is only C 1,1 . We can further prove that the M-D transonic shock-front solution is 
stable with respect to the C 2 ' a supersonic perturbation. 

When the regularity of the steady perturbation is C 3 ' a or higher : 1 1 ip ~ — ip^ \ g d ~ ^ < 
a, we have introduced another simpler approach, the implicit function approach, to 
deal with the existence and stability problem in Chen-Feldman [47] . 

Other approaches can be found in Canic-Keyfitz-Lieberman [31], Canic-Keyfitz- 
Kim [33], Chen-Chen- Feldman [43], Chen [72, 73], Zheng [215], and the references 
cited therein. 
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5.5.2. Weak Convergence Approaches. Recent investigations on conservation laws 
based on weak convergence methods suggested that the method of compensated 
compactness be amenable to flows which exhibit both elliptic and hyperbolic regimes. 
In [163] (also see [164]), Morawetz layed out a program for proving the existence 
of the steady transonic flow problem about a bump profile in the upper half plane 
(which is equivalent to a symmetric profile in the whole plane). It was shown that, 
if the key hypotheses of the method of compensated compactness could be satis- 
fied, now known as a "compactness framework" (see Chen [39]), then indeed there 
would exist a weak solution to the problem of flow over a bump which is exhibited 
by subsonic and supersonic regimes, i.e., transonic flow. 

The compactness framework for system (5.17) can be formulated as follows. Let 
a sequence of functions w e (x) = (u e ,v £ )(x) defined on an open set f2 C M 2 satisfy 
the following set of conditions: 

(A.l) q £ (x) = |w £ (x)| < <7(j a.e. in £1 for some positive constant q$ < q ca v', 

(A. 2) V x • Q±(w e ) are confined in a compact set in H^(Q) for entropy pairs 
Q± = (Qi±iQ2±), where Q±(w £ ) arc confined to a bounded set uniformly in 

When (A.l) and (A. 2) hold, the Young measure f x (w),w = (u,v), determined 
by the uniformly bounded sequence of functions w e is constrained by the following 
commutator relation: 

(fx, Ql + Q2- -Q1-Q2+) = (fx, <9l + )(f X , Q2-) - (fx, Ql-)(fx, Q2+)- (5.21) 

The main point for the compensated compactness framework is to prove that f x is a 
Dirac mass by using entropy pairs, which implies the compactness of the sequence 
w e (x) = (u e ,w e )(x) in L 1 1 oc (f2). In this context, It is needed Morawetz [163] to 
presume the existence of an approximating sequence parameterized by e to their 
problems satisfying (A.l) and (A. 2) so that they could exploit the commutator 
identity and obtain the strong convergence in L loc (Q) to a weak solution of their 
problems. 

As it turns out, there is a classical problem where (A.l) and (A. 2) hold trivially, 
i.e., the sonic limit of subsonic flows. In that case, we return to the result by Bers 
[15] and Shiffman [189], which says that, if the speed at infinity, goo, is less than some 
q, there is a smooth unique solution of the problem and ask what happens as q^ /* 
q. In this case, the flow develops sonic points and the governing equations become 
degenerate elliptic. Thus, if we set e = q — qca and examine a sequence of exact 
smooth solutions to our system, we see trivially that (A.l) is satisfied since \q £ \ < g*, 
and (A. 2) is also satisfied since V x • Q±(w £ ) = along our solution sequence. The 
effort is in finding entropy pairs which guarantee the Young measure ^ x reduces to 
a Dirac mass. Ironically, the original conservation equations of momentum in fact 
provide two sets of entropy pairs, while the irrotationality and mass conservation 
equations provide another two sets. This observation has been explored in detail 
in Chen-Dafermos-Slemrod-Wang [45]. 

What then about the fully transonic problem of flow past an obstacle or bump 
where qoo > q! In Chen-Slemrod-Wang [61], we have provided some of the in- 
gredients to satisfying (A.l) and (A. 2) as explained in §5.4. On the other hand, 
(A. 2) is easily obtained from the viscous formulation by using a special entropy 
pair of Osher-Hafez- Whitlow [172]. In fact, this entropy pair is very important: 
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It guarantees that the inviscid limit of the above viscous system satisfies a physi- 
cally meaningful "entropy" condition (Theorem 2 of [172]). With (A.l) and (A. 2) 
satisfied, the compensated compactness argument then applies, to yield the strong 
convergence in Lj oc (Q) of our approximate solutions. 

Some further compensated compactness frameworks have been developed for 
solving the weak continuity of solutions to the Gauss-Codazzi-Ricci equations in 
Chen-Slemrod-Wang [62, 63]. Also see Dacorogna [88], Evans [103], and the refer- 
ences cited therein for various weak convergence methods and techniques. 

6. Shock Reflection-Diffraction and Self-Similar Solutions 

One of the most challenging problems is to study solutions with data that give 
rise to self-similar solutions (such solutions include Riemann solutions). For the 
Euler equations (1.4) for x <G R 2 , the self-similar solutions 

(P,v,p) = (p,v,p)(£,r)), {t,v)=x/t, (6-1) 

are determined by 

d s ( P u) + d v ( P v) = -2 P , 

ddpU 2 + p) + d n ( P UV) = -3 P U, 

dz(puv) + d v (pv 2 +p) = -3pv 1 

ds(pU(E + p/p)) + d n { P V{E + p/p)) = -2p(E + p/p), 

where (£/, V") = v - (£, 77) is the pseudo-velocity and E = e + \{U 2 + V 2 ). 
The four eigenvalues are 

. V UV±cVU 2 + V 2 -c 2 

\ = — (2-multiphcity), Aj 



v v - •>» - U 2 -c 2 

where c = \Jp p {p, S) is the sonic speed. When U 2 + V 2 > c 2 , system (6.2) is 
hyperbolic with four real eigenvalues and the flow is called pseudo-supersonic, sim- 
ply called supersonic without confusion. When U 2 + V 2 < c 2 , system (6.2) is 
hyperbolic-elliptic composite type (two repeated eigenvalues are real and the other 
two are complex): two equations are hyperbolic and the other two are elliptic. The 
region U 2 + V 2 = c 2 in the (£, ?y)-plane is simply called the sonic region in the 
flow. In general, system (6.2) is hyperbolic-elliptic mixed and composite type, and 
the flow is transonic. For a bounded solution (p,v,p), the flow must be supersonic 
when £ 2 + rj 2 — > 00. 

An important prototype problem for both practical applications and the theory 
of M-D complex wave patterns is the problem of diffraction of a shock-front which 
is incident along an inclined ramp. When a plane shock-front hits a wedge head 
on, a self-similar reflected shock-front moves outward as the original shock-front 
moves forward (cf. [13, 34, 84, 116, 165, 183, 201]). Then the problem of shock 
reflection-diffraction by a wedge can be formulated as follows: 

Problem 6.1 (Initial-boundary value problem). Seek a solution of system 
(1.4) satisfying the initial condition at t = 0: 

v _ j(po,0,0,po), \xi\ > x\t<m6 w ,xi > 0, 
\(Po,ui,0,pi), xi < 0; 

and the slip boundary condition along the wedge boundary: 

v • n = 0, (6.4) 
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where n is the exterior unit normal to the wedge boundary, and state (0) and (1) 
satisfy 



(Pi - Po){pi - Po) Pi (7 + l)pi - (7- l)Po n . n (Rt ., 
V PoPi Po (7 + l)Po - (7- IjPi 

Given po,po,pi, and 7 > 1, the other variables u\ and p\ are determined by 
(6.5). In particular, the Mach number M\ = ^ for state (1) is determined by 

A/r2 _ 2(pi-p ) 2 
JU 1 " po((7+l)Pi-(7-l)Po)' 



X2 




Figure 1 . Initial-boundary value problem 



Since the initial-boundary value problem, Problem 6.1, is invariant under the self- 
similar scaling, we seek self-similar solutions (6.1) governed by system (6.2). Since 
the problem is symmetric with respect to the axis 77 = 0, it suffices to consider the 
problem in the half-plane 77 > outside the half- wedge: 

A := {£ < 0,7] > 0}U {77 > £t&n6 w , £ > 0}. 

Then Problem 6.1 in the (t, x)~coordinates can be formulated as the following 
boundary value problem in the self-similar coordinates (£,?/): 

Problem 6.2 (Boundary value problem in the unbounded domain). 

Seek a solution to system (6.2) satisfying the slip boundary condition on the wedge 
boundary: (U, V) ■ n = on dk = {£ < 0, 77 = 0} U {£ > 0, n > £tan6> w }, the 
asymptotic boundary condition as £ 2 + rj 2 — > 00: 



{p,U + i,V + n,p) 



(Po,0,0,p ), C > ^o,?7 > Ctan^ 
(pi,ui,0,pi), C < ^o, ?7 > 0. 



For our problem, since <pi does not satisfy the slip boundary condition (6.11), 
the solution must differ from ip\ in {£ < ^o}H A, thus a shock diffraction-diffraction 
by the wedge vertex occurs. The experimental, computational, and asymptotic 
analysis shows that various patterns of reflected shock-fronts may occur, including 
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Figure 2. Boundary value problem in the unbounded domain A 

the regular and Mach reflections (cf. [13, 84, 116, 121, 128, 129, 132, 133, 148, 155, 
165, 191, 192, 197, 198, 199, 201, 209]). It is expected that the solutions of Problem 
6.2 contain all possible patterns of shock reflection-diffraction configurations as 
observed. For a wedge angle 8 W € (0, ■§), different reflection-diffraction patterns 
may occur. Various criteria and conjectures have been proposed for the existence 
of configurations for the patterns (cf. Ben-Dor [13]). One of the most important 
conjectures made by von Neumann [170, 171] in 1943 is the detachment conjecture, 
which states that the regular reflection-diffraction configuration may exist globally 
whenever the two-shock configuration (one is the incident shock-front and the other 
the reflected shock- front) exists locally around the point Pq (see Fig. 3). The 
following theorem was rigorously shown in Chang- Chen [34] (also see Sheng-Yin 
[187], Bleakney-Taub [20], Neumann [170, 171]). 

Theorem 6.1 (Local theory). There exists 9d = &d(Po, Pi,l) € (0, § ) such that, 
when 9 W € (0d, §), there are two states (2):(p 2 , U 2 , V 2 ,p 2 ) an d (p!> ^2> suc ^ 
that \{U$,V 2 a )\ > \{U b 2 , V 2 b )\ and \(U b , V 2 b )\ < c{p b , S b 2 ). 

The von Neumann Detachment Conjecture ([170, 171]): There exists a 
global regular reflection- diffraction configuration whenever the wedge angle 9 W is in 

It is clear that the regular reflection-diffraction configuration is not possible 
without a local two-shock configuration at the reflection point on the wedge, so this 
is the weakest possible criterion. In this case, the local theory indicates that there 
are two possible choices for state (2). There had been a long debate to determine 
which one is more physical for the local theory; see Courant-Friedrichs [84], Ben- 
Dor [13], and the references cited therein. Since the reflection-diffraction problem 
is not a local problem, we take a different point of view that the selection of state 
(2) should be determined by the global features of the problem, more precisely, by 
the stability of the configuration with respect to the wedge angle 6 W , rather than 
the local features of the problem. 

Stability Criterion to Select the Correct State (2) (Chen-Feldman [49]): 

Since the solution is unique when the wedge angle 9 W = ^, it is required that the 
global regular reflection- diffraction configuration be stable and converge to the unique 
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normal reflection solution when 6 W — > 5, provided that such a global configuration 
can be constructed. 

We employ this stability criterion to conclude that our choice for state (2) must 
be {p2iU%,V£,p%). In general, (p|, U$ , V^pf) mav be supersonic or subsonic. If 
it is supersonic, the propagation speeds are finite and state (2) is completely deter- 
mined by the local information: state (1), state (0), and the location of the point 
Pq. This is, any information from the reflected region, especially the disturbance 
at the corner P 3 , cannot travel towards the reflection point Pq. However, if it is 
subsonic, the information can reach Pq and interact with it, potentially altering the 
reflection-diffraction type. This argument motivated the following conjecture: 

The von Neumann Sonic Conjecture [170, 171]: There exists a regular 
reflection- diffraction configuration when 9 W € (# s , § ) for 9 S > Qd such that | (C/f , ^2°)! > 
cl at P . 




Figure 3. Regular reflection-diffraction configuration 

If state (2) is sonic when 9 W = 6 S , then \(U§,V£)\ > c% for any 9 W <E (9 S , § ). 
This sonic conjecture is stronger than the detachment one. In fact, the regime 
between the angles 9 S and 9^ is very narrow and is only fraction of a degree apart; 
see Sheng-Yin [187]. 

Following the argument in Chen-Feldman [50], we have 

Theorem 6.2 (Chen-Feldman [50]). Let (p,U,V,p) be a solution of Problem 6.2 
such that (p,U,V,p) is C 0,1 in the open region P0P1P2P3 and the gradient of the 
tangential component of (U, V) is continuous across the sonic arc T SO nic- Let £l\ be 
the subregion of Q formed by the fluid trajectories past the sonic arc T son i C . Then, 
in ili, the potential flow equation for self-similar solutions: 

V- (p(Vtp,<p)Vtp) +2p(V(^) = 0, (6.6) 

with pdVt^l 2 , 99) = (pq 1 — (7 — l)(<p + 5IV93I 2 )) 0-1 , coincides with the full Euler 
equations (6.2), that is, equation (6.6) is exact in the domain fii. 

The regions such as f?i also exist in various Mach reflection-diffraction config- 
urations. Theorem 6.2 applies to such regions whenever the solution (p,U,V,p) is 
C 0,1 and the gradient of the tangential component of (U, V) is continuous. In fact, 
Theorem 6.2 indicates that, for the solutions <p of (6.6), the C^-regularity of <p 
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and the continuity of the tangential component of the velocity field (U, V) = Vip 
are optimal across the sonic arc r so „i C . 

Equation (6.6) is a nonlinear equation of mixed elliptic-hyperbolic type. It is 
elliptic if and only if \Vip\ < c(\W(p\ 2 , (p, pi ), which is equivalent to 



|Vp| << Z 4^po,7):=y^- T (pr 1 -(7-l)^)- (6-7) 

For the potential equation (6.6), shock-fronts are discontinuities in the pseudo- 
velocity Vip. That is, if D + and D~ := D \ D + are two nonempty open subsets 
of D C K 2 , and S :— dD + n D is a C 1 -curve where Dip has a jump, then ip e 
w iol ( D ) n Cl ( D± U 5) n C 2 (L> ± ) is a global weak solution of (6.6) in D if and only 
if tp is in W l( £°(D) and satisfies equation (6.6) in D ± and the Rankinc-Hugoniot 
condition on S: 

[ P (|V^| 2 ,^)V^-n] 5 =0. (6.8) 
Then the plane incident shock-front solution in the (t, x)-coordinates with states 
(p, V x 3>) = (po, 0, 0) and (pi, iti, 0) corresponds to a continuous weak solution ip of 
(6.6) in the self-similar coordinates (C, 77) with the following form: 

MZ,V) = -\(Z 2 +V 2 ) for C>Co, (6.9) 

Vi(t,ri) = -\(e+V 2 )+Mt-to) for C<Co, (6.10) 



respectively, where «i = y 2(p \^{l Pl ] rp $ ^ > and Co = > are the ve- 

locity of state (1) and the location of the incident shock- front, uniquely determined 
by (p ,pi,j) through (6.8). Then P = (Co, Co tan&u,) in Fig. 2, and Problem 6.2 
in the context of the potential flow equation can be formulated as 

Problem 6.3 (Boundary value problem) (cf. Fig. 2). Seek a solution <p of (6.6) 
in the self-similar domain A with the boundary condition on dA: 

V(^-n| 9A = 0, (6.11) 

and the asymptotic boundary condition at infinity: 

jVo for C > Co, V > C tan W , 2 2 

, t ^ t n when +r] 2 ^ 00, (6.12) 

[<^i /or C < Co, 1? > 0, 

where (6.12) holds in the sense that lim \\ip — ^||c(A\s R (o)) = 0- 

In Chen-Feldman [49], we first followed the von Neumann criterion and the 
stability criterion introduced above to establish a local existence theory of regular 
shock reflection near the reflection point Po m the level of potential flow, when the 
wedge angle is large and close to ^. In this case, the vertical line is the incident 
shock- front that hits the wedge at P = (Co, Co tan 6 W ), and state (0) and state 
(1) ahead of and behind S are given by ipo and ipi defined in (6.9) and (6.10), 
respectively. The solutions ip and ipi differ only in the domain P0P1P2P3 because 
of shock diffraction by the wedge, where the curve P0P1P2 is the reflected shock- 
front with the straight segment PoPl. State (2) behind P0P1 can be computed 
explicitly with the form: 

¥> 2 (C, V) = -^(C 2 + V 2 ) + « 2 (C - Co) + (V - Co tan tan^, (6.13) 
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which satisfies V<p • n = on <9A n {£ > 0}; the constant velocity 112 and the angle 
6 S between PqPi and the £-axis are determined by (0 W , po, pi, 7) from the two 
algebraic equations expressing (6.8) and continuous matching of states (1) and (2) 
across P$P\ as in Theorem 6.1. Moreover, p 2 is the unique solution in the domain 
P0P1P4, as argued in [34, 183]. Hence P1P4 := T sonic = dfl fl dB C2 (u 2 ,U2ta,n9 w ) 
is the sonic arc of state (2) with center (w 2 , 112 tan^) and radius c 2 . 

It should be noted that, in order that the solution ip in the domain f2 is a 
part of the global solution to Problem 6.3, that is, if satisfies the equation in the 
sense of distributions in A, especially across the sonic arc P1P4, it is required that 
V(<^ — 1P2) ■ n |piP 4 = 0. That is, we have to match our solution with state (2), 
which is the necessary condition for our solution in the domain Q to be a part of 
the global solution. To achieve this, we have to show that our solution is at least 
C 1 with V(<p — (P2) = across P\Pa- Then Problem 6.3 can be reformulated as the 
following free boundary problem: 

Problem 6.4 (Free boundary problem). Seek a solution <p and a free bound- 
ary P1P2 = {£ = f(T])} suc -h that 

(i) / G C l ' a and 

^={?>/W}nA={^< W }nA; (6.14) 

(ii) ip satisfies the free boundary condition (6.8) along P1P2; 

(iii) (p G C 1,a (fl) n C 2 (fl) solves (6.6) in fl, is subsonic in fl, and satisfies 

(ip-ip2,V{ip-ip2)-n)\ Pl p 2 =0, (6.15) 
Vp • n|p 3 p 4 ur syrom = 0. (6.16) 

The boundary condition on T symm implies that /'(0) = and thus ensures the 
orthogonality of the free boundary with the £-axis. Formulation (6.14) implies 
that the free boundary is determined by the level set ip = <p\. The free boundary 
condition (6.8) along P1P2 is the conormal boundary condition on P\Pi- Condition 
(6.15) ensures that the solution of the free boundary problem in SI is a part of 
the global solution. Condition (6.16) is the slip boundary condition. Problem 6.4 
involves two types of transonic flow: one is a continuous transition through the 
sonic arc P1P4 as a fixed boundary from the supersonic region (2) to the subsonic 
region fl; the other is a jump transition through the transonic shock-front as a free 
boundary from the supersonic region (1) to the subsonic region Q. 

In Chen-Feldman [49, 50, 51], we have developed a rigorous mathematical ap- 
proach to solve the von Neumann sonic conjecture, Problem 6.4, and established 
a global theory for solutions of regular reflection-diffraction up to the sonic angle, 
which converge to the unique solution of the normal shock reflection when 9 W tends 
to 7r/2. For more details, see Chen-Feldman [49, 50, 51]. 

The mathematical existence of Mach reflection-diffraction configurations is still 
open; see [13, 74, 84]. Some progress has been made in the recent years in the 
study of the 2-D Riemann problem for hyperbolic conservation laws; see [35, 36, 
39, 114, 115, 116, 122, 134, 140, 143, 144, 182, 184, 186, 188, 194, 213, 214] and the 
references cited therein. 

Some recent developments in the study of the local nonlinear stability of mul- 
tidimensional compressible vortex sheets can be found in Coulombel-Secchi [83], 
Bolkhin-Trakhinin [20], Trakhinin [200], Chen- Wang [68], and the references cited 
therein. Also see Artola-Majda [7]. For the construction of the non-self-similar 
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global solutions for some multidimensional systems, see Chen- Wang- Yang [67] and 
the references cited therein. 

7. Divergence-Measure Vector Fields and Multidimensional 

Conservation Laws 

Naturally, we want to approach the questions of existence, stability, uniqueness, 
and long-time behavior of entropy solutions for M-D hyperbolic conservation laws 
with neither specific reference to any particular method for constructing the so- 
lutions nor additional regularity assumptions. Some recent efforts have been in 
developing a theory of divergence-measure fields to construct a global framework 
for the analysis of solutions of M-D hyperbolic systems of conservation laws. 

Consider system (1.1) in R d . As discussed in §2.5, the BV bound generically 
fails for the M-D case. In general, for M-D conservation laws, especially the Eulcr 
equations, solutions of (1.1) are expected to be in the following class of entropy 
solutions: 

(i) u(t,x) e M(R d + +1 ), or LP(R d + +1 ), 1 < p < oo; 

(ii) u(t, x) satisfies the Lax entropy inequality: 

fir, := d tV (u(t, x)) + V x • q(u(t, x)) < (7.1) 

in the sense of distributions for any convex entropy pair (77, q) : K™ — > KxIR d 
so that r)(u(t, x)) and q(u(t, x)) are distributions. 

The Schwartz Lemma infers from (7.1) that the distribution fi v is in fact a Radon 
measure: div(i >x ) (77(11^, x)),q(u(f, x))) e A4(R+ +1 ). Furthermore, when u e L°°, 
this is also true for any C 2 entropy-entropy flux pair (77, q) (77 not necessarily convex) 
if (1.1) has a strictly convex entropy, which was first observed in Chen [38]. More 
generally, we have 

Definition. Let V C M. N be open. For 1 < p < 00, F is called a X>M p (X>)-field if 
F e LP(V;R N ) and 

\\F\\vmp(v) : = \\ f \\lp(v-m n ) + \\^ vF \\m(v) < °°; (7-2) 

and the field F is called a VM ext (D)-field if F G M(D;R N ) and 

||F|b M -*(x>) := ||(F,divF)|| A4 (x)) < 00. (7.3) 

Furthermore, for any bounded open set V C K w , F is called a 2?Mf oc (R A ')-field if 
F e VM P (V); and F is called a VM^(R N )-fteld if F e VM ext (V). A field F is 
simply called a X>M-field in V if F e VM P {V), 1 < p < 00, or F e VM ext (V). 

It is easy to check that these spaces, under the respective norms ||F||d MP (x>) 
and ||F||x)M' !:ct (X')! are Banach spaces. These spaces are larger than the space of 
BV-fields. The establishment of the Gauss-Green theorem, traces, and other prop- 
erties of -BV functions in the 1950s (cf. Federer [105]; also Ambrosio-Fusco-Pallara 
[5], Giusti [112], and Volpert [203]) has significantly advanced our understanding 
of solutions of nonlinear PDEs and related problems in the calculus of variations, 
differential geometry, and other areas, especially for the 1-D theory of hyperbolic 
conservation laws. A natural question is whether the 2?M-fields have similar prop- 
erties, especially the normal traces and the Gauss-Green formula to deal with en- 
tropy solutions for M-D conservation laws. 
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On the other hand, motivated by various nonlinear problems from conservation 
laws, as well as for rigorous derivation of systems of balance laws with measure 
source terms from the physical principle of balance law and the recovery of Cauchy 
entropy flux through the Lax entropy inequality for entropy solutions of hyperbolic 
conservation laws by capturing the entropy dissipation, a suitable notion of normal 
traces and corresponding Gauss-Green formula for divergence-measure fields are 
required. 

Some earlier efforts were made on generalizing the Gauss-Green theorem for 
some special situations, and relevant results can be found in Anzcllotti [6] for an 
abstract formulation for F e L°° , Rodrigues [178] for F e L 2 , and Ziemer [216] 
for a related problem for divF <G L 1 ; also see Baiocchi-Capelo [10] and Brezzi- 
Fortin [29]. In Chen-Frid [54], an explicit way to calculate the suitable normal 
traces was first observed for F G DM 00 , under which a generalized Gauss-Green 
theorem was shown to hold, which has motivated the development of a theory of 
divergence-measure fields in [54, 64, 66]. 

Some entropy methods based on the theory of divergence-measure fields pre- 
sented above have been developed and applied for solving various nonlinear prob- 
lems for hyperbolic conservation laws and related nonlinear PDEs. These problems 
especially include (i) Stability of Riemann solutions, which may contain rarefaction 
waves, contact discontinuities, and/or vacuum states, in the class of entropy solu- 
tions of the Euler equations in [42, 54, 55]; (ii) Decay of periodic entropy solutions 
in [52]; (hi) Initial and boundary layer problems in [54, 60, 64, 202]; (iv) Rigorous 
derivation of systems of balance laws from the physical principle of balance law and 
the recovery of Cauchy entropy flux through the Lax entropy inequality for entropy 
solutions of hyperbolic conservation laws by capturing the entropy dissipation in 
[66]; 

It would be interesting to develop further the theory of divergence-measure fields 
and more efficient entropy methods for solving more various problems in PDEs and 
related areas whose solutions are only measures or LP functions. For more details, 
see [54, 64, 66]. 

We also refer the reader to the related papers on other aspects of multidimen- 
sional conservation laws in this volume. 
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